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Abstract 

This document is a companion for the Maple program Discrete 
series and T^T-types for U(p, q) available on 

\protect\vrule widthOpt\protect\href {http : //www. math. jussieu. fr/\string~vergne/}{] 

We explain an algorithm to compute the multiplicities of an irre- 
ducible representation of U(p) x U(q) in a discrete series of U(p,q). 
It is based on Blattner's formula. We recall the general mathematical 
background to compute Kostant partition functions via multidimen- 
sional residues, and we outline our algorithm. We also point out some 
properties of the piecewise polynomial functions describing multiplic- 
ities based on Paradan's results. 
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Introduction 

The present article is a user's guide for the Maple program Discrete series and 
K types for U(p,q), available at 

\protect\vrule widthOpt\protect\href {http : //www .math. jussieu . f r/\string~vergne/}{http : , 

In the first part, we explain what our program does with simple examples. 
The second part sets the general mathematical background. We recall Blattner's 
formula, and we discuss the piecewise (quasi)-polynomial behavior of multiplicities 
for a discrete series of a reductive real group. In the third part, we outline the 
algorithm of computing partition functions for arbitrary set of vectors, based on 
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Jeffrey-Kirwan residue and maximally nested subsets. In the fourth part, we spe- 
cialize the method to U (p, q) . We in particular explain how to compute maximally 
nested subsets of the set of non compact positive roots of U (p, q) . 

In the last two sections, we give more examples and some details on the im- 
plementation of the algorithms for our present application. 

Here G is the Lie group G = U(p, q) and K = U(p) X U(q) is a maximal 
compact subgroup of G. Of course all these issues can be addressed for the other 
reductive real Lie groups following the same approach described here. To introduce 
the function we want to study, we briefly establish some notations (see Sec. 12. 2\i . 
We denote by tt x a discrete series representation with Harish Chandra parameter 
A. The restriction of 7r A to the maximal compact subgroup K decomposes in irre- 
ducible finite dimensional K representations with finite multiplicities, in formula 

^ A \k = Yl m W 

where we sum over K: the classes of irreducible finite dimensional K represen- 
tations, [i being the Harish Chandra parameter of r^. 

m x is a finite number called the multiplicity of the K-type r^, or simply of /x, in 
7r A and this paper addresses the question of computing m^. 

The algorithm described in this paper checks whether a certain .fT-type r M appears 
in the if-spectrum of a discrete series ir x by computing the multiplicity of such a 
-ftT-type. The input /i can also be a symbolic variable as we will explain shortly. 
By Blattner's formula, computing the K-multiplicity is equivalent to compute the 
number of integral points, that is a partition function, for specific polytopes. We 
thus use the formulae developed in [3] to write the algorithm. 

One important aspect of these results is that the input datas (A, /i) can also be 
treated as parameters. Thus, in principle, given (Xq,/j,q), we can output a convex 
cone, containing (Ao, Ho), and a polynomial function of the parameters (A, n) whose 
value at (A, /i) is the multiplicity of the K-type /x in the representation 7r A as long as 
we stay within the region described by the cone and (A, //) satisfy some integrality 
conditions to be defined later. We also plan to explicitly decompose our parameter 
space (A, /j) in such regions of polynomiality in a future study, thus describing fully 
the piecewise polynomial function at least for some low rank cases. 

In practise here, we only address a simpler question: we will fix A, /i and 
a direction v and compute a piecewise polynomial function of t coinciding with 
m ^,+tv on integers t. This way, we can also check if a direction v is an asymptotic 
direction of the i"C-spectrum, in the sense explained in Sec. 13.7.31 

The Atlas of Lie Groups and Representations, [I], within the problem of clas- 
sifying all of the irreducible unitary representations of a given reductive Lie group, 
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addresses in particular the problem of computing K-types of discrete series. The 
multiplicities results needed for the general unitary problem is of different nature 
as we are going to explain. 

Given as input A and some height h, Atlas computes the list with multiplic- 
ities of all the representations occurring in n x of height smaller than h. But the 
efficiency, in this setting, is limited by the height. In contrast, the efficiency of our 
program is unsensitive to the height of A, fj,, but the output is one number: the 
multiplicity of [i in 7r A . It takes (almost) the same time to compute the multiplicity 
of the lowest if -type of tt x (fortunately the answer is 1) than the multiplicity of a 
representation of very large height. Our calculation are also very sensitive to the 
rank: p + q — 1. 

For other applications (weight multiplicities, tensor products multiplicities) 
based on computations of Kostant partition functions in the context of finite di- 
mensional representations, see [4], jfj], [2]. 

1 The algorithm for Blattner's formula: main 
commands and simple examples 

Let p,q be integers. We consider the group G = U(p,q). The maximal compact 
subgroup is K := U(p) x U(q). More details in parametrization are given in Section 

m 

A discrete series representation 7r A is parametrized according to Harish-Chandra 
parameter A, that we input as discrete: 

discrete := [[Ai, . . . , X p ], [71 , . . . , 7,]] . 

Here A,, jj are integers if p + q is odd, or half-integers if p + q is even. They 
are all distinct. Furthermore Ai > • • • > A p and 71 > ■ • • > 7 g . 

A unitary irreducible representation of K (that is a couple of unitary irreducible 
representations of U(p) and of U(q)) is parametrized by its Harish-Chandra pa- 
rameters that we input as Krep: 

Krep := [[ai,a 2 , ■ ■ ■ ,a p ], [b l7 . . .,b q ]] 

with a\ > ■ ■ ■ > a p and b\ > ■ ■ ■ > b q . 

Here are integers if p is odd, half-integers if p is even. Similarly bj are 
integers if q is odd, half-integers if q is even. 

As we said, our objective is to study the function for \x € K where /j, = Krep 
and A = discrete. 
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The examples are runned on a MacBook Pro, Intel Core 2 Duo, with a Processor 
Speed of 2.4 GHz. The time of the examples is computed in seconds. They are 
recorded by 

TT 

Some of these examples are very simple and can be checked by hand (as we 
did, to reassure ourselves). Other examples are given at the end of this article. 

To compute the multiplicity of the if -type given by Krep in the discrete series 
with parameter given by discrete, the command is 

>discretemult (Krep , discrete , p , q) 
Example 1 

Krep: =[[207/2, -3/2], [3/2, -207/2]]; 
discrete : = [ [5/2, -3/2], [3/2, -5/2]]; 

>discretemult (Krep, discrete ,2,2); 

101 

Here is another example of that our program can compute. 
Example 2 

We consider the discrete series indexed by 

lambda33:=[[ll/2, 7/2, 3/2], [9/2, 5/2, 1/2]]; 
Its lowest K-type is 

lowlambda33 : = [[7,4,1], [5,2,-1]]; 

Of course, the multiplicity of the lowest K-type is 1. Our programm fortunately returns the value 1 in 
0.03 seconds. 

Consider now the representation of K with parameter: 

biglambda33 : = [[10006, 4, -9998], [10004, 2, -10000]]; 
Then the multiplicity of this ii"-type is computed computed in 0.05 seconds as 
> discretemult(biglambda33,lambda33,3,3) ; 
2500999925005000; 
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Here are two other examples that verify the known behavior of holomorphic dis- 
crete series. The notation ab... that we use to label the discrete series parameters 
is introduced in 14.11 and it very effective to picture the situation, but it is not 
relevant to understand the following computation. 

Example 3 

Consider a holomorphic discrete series of type "aaabbb" for G = U(3, 3) (see Sec l4. Il l with lowest i\-type 
of dimension 1 . We verify that the multiplicity is 1 for fj, in the cone spanned by strongly orthogonal non 
compact positive roots. 

hol33:=[ [11/2, 9/2, 7/2] , [5/2,3/2,1/2]] ; 
lowhol33:=[[7, 6, 5], [1, 0, -1]]; 

bighol33:=[[7+1000, 6+100, 5+10], [1-10, -100, -1-1000]]; 

>discretemult (lowhol33 , hol33 ,3,3); 

1 

TT:= 0.052 

>discretemult (bighol33 , hol33 , 3 , 3) ; 
1 

TT:= 1.003 

Example 4 

Consider a holomorphic discrete series of type "aaabbb" for G = U(3, 3) (see Sec. I4.1[l with lowest _R"-type 
of dimension d. We verify that the multiplicities are bounded by d. 

Hol33 : = [ [27/2 ,9/2, 7/2], [5/2, 3/2, -5/2]]; 
lowHol33: = [[15, 6, 5], [1, 0, -4]]; 

bigHol33:=[ [15+1000, 6+1000, 5+1000], [1-1000, -1000, -4-1000]]; 
verybigHol33:=[ [15+100000, 6+10000, 5+10000], [1-10000, -10000, -4-100000]]; 

>discretemult (lowHol33 ,Hol33 ,3,3); 

1 

TT:= 0.069 

>discretemult (bigHol33 ,Hol33 ,3,3); 
4 

TT: =0.971 

>discretemult ( verybigHol33 , Hol33 , 3 , 3) ; 
4 

TT:= 0.873 



Fix now a K-type //, a direction v given by a dominant weight for K ( more 
details in Sec. I3.T.3H and a discrete series parameter A. The half-line [i + tv stays 
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inside the dominant chamber for K. A very natural question is that of investigating 
the behavior of the multiplicity function as a function of t £ Z, when we move 
from \i along the positive v direction, that is the function t — > m^ +t ^, t > 0. 
The answer to this question will be given by two sets of datas: a covering of N 
determined by a finite number of closed intervals ij C K with integral end points, 
that is N = Ui<j< s (/j nN), together with polynomial functions Pj(i), 1 < i < s, of 
degree bounded by pq — (p + q — 1) , that compute the multiplicity on such intervals 
Ii\ in formula m^ +t ^ = Pi(t) for t 6 ij fl N. 

We remark two aspects. First the {ij CI N, 1 < i < s} constitutes a covering 
in the sense that we recover all of N but (Ij n N) n (ij+i H N) can intersect in 
the extreme points and hence in this case Pj and Pj+i have to coincide on the 
intersection. Secondly the "intervals" Ii HN can be reduced just to a point, so that 
the polynomial Pi (if not constant) is not uniquely determined by its value on one 
point !. More generally, if the length of the interval Ii is smaller that the degree 
of Pj, the polynomial Pi is not uniquely determined. 

Because of our focus on the polynomiality aspects we keep in the output of the 
algorithm the polynomials Pj even if the intervals Ii are reduced to a point (or 
with small numbers of integral points). 

In our application, fi will be the lowest K-type and we will give some examples of 
the situations occurring, in particular we will examine the following cases: 

• The first example outputs a covering of N given by a unique interval and 
a polynomial function on N that computes the multiplicity. Thus in this 
case, all the integer points on the half line [i + tv, t > give rise to K-types 
that appear in the restriction of the discrete series, in particular v is an 
asymptotic direction, (see Sec. 12.41 and |3.T.3[) . 

• In the other examples, the covering of N has at least two intervals and 
illustrate different situations. 

To compute, in the sense we just explained, the multiplicity of the K-type fi + tv 
in the discrete series with parameter discrete, [i being the lowest if -type moving 
in the positive direction v, labeled by direction , the command is 

>f unction_discrete_mul_direction_lowest (discrete .direct ion, p,q) ; 



Example 5 
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discrete:=[[5/2, -3/2], [3/2, -5/2]]; 
direction:=[[l,0] , [0,-1]] ; 

>function_discrete_mul_direction_lowest (discrete , [ [1 ,0] , [0 , -1] ] , 2,2); 
[[t+1, [0, inf]]] 

(inf stands for oo ). 

Here the covering of N is N = [0, oo] n N and the polynomial is P(t) =t + l. The 
output explicitly compute: 

m£ + ^ = i + l, teN, t>0 

with \i = [[7/2, —3/2], [3/2, —7/2]] the lowest if -type of the representation 7r A (See 
Ex H] for computation of the lowest if- type). Thus we compute the multiplicity, 
starting from the lowest if-type when we are moving off it in the direction of v. 
In particular for t = 100 we get "i^ +10 (h; = 101 as predicted in ExfTJ since 

discrete:=[[5/2, -3/2], [3/2, -5/2]] 

and // + WOv is equal to 

Krep: = [ [207/2, -3/2], [3/2, -207/2]]; 

Example 6 

discrete: = [[9, 7], [-1, -2, -13]]; 
directions [[1, 0], [0, 0, -1]]; 

>f unction_discrete_mul_direction_lowest (discrete , direction, 2,3) ; 
[[l+(l/2)*t-(l/2)*t"2, [0, 0]], [1, [1, inf]]] 

Here the covering is N = ([0, 0] n N) U ([1, oo] n N and the polynomials are P\(t) = 
1 + \t - \t 2 on [0, 0] n N and P 2 (t) = 1 on ([1, oo] n N. Observe that [0, 0] n N = [0] 
is just a point and that -Pi(O) = 1 as it should be, since \x is the lowest K-type. 
Explicitly we simply compute 

m^=l, t€N, t>0. 

We conclude with one example in which the multiplicity grows: the last poly- 
nomial is not zero and has degree two. We give more examples at the end of this 
article. 
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Example 7 



discrete:= [[57/2, 39/2, 3/2], [51/2, 5/2, -155/2]]; 
direction: = [[1, 0, 0], [0, 0, -1]]; 

>function_discrete_mul_direction_lowest (discrete, direct ion, 3, 3) ; 

[(l/24)*t-4+(5/12)*t-3+(35/24)*t-2+(25/12)*t+l, [0, 16]], [-3059+ (2242/3) *t-( 133/2) *t'2+(19/6) *t'3 , [17, 21]], 
[-11914+(9597/4)*t-(4367/24)*t-2+(27/4)*t-3-(l/24)*f4, [22, 40]], [100016-8664*t+228*t"2 , [41, inf ] ] 



That is for A =discrete and [i = [[30, 20, 1], [26, 2, -79]] the lowest K-type 



' (1/24) *t 4 + (5/12) *t 3 + (35/24) *t 2 + (25/12) *t + l < t < 16 

-3059 + (2242/3) * t - (133/2) * t 2 + (19/6) *t 3 17 < t < 21 

-11914 + (9597/4) * t - (4367/24) * t 2 + (27/4) * t 3 - (1/24) * t 4 22 < t < 40 

100016 - 8664 * t + 228 * t 2 t > 41 



The time to compute the example is TT := 0.835 and the formula says for instance 
that mj +2000000 - = 91198 2 6 72 1000 16 

To compare with other parametrizations of the discrete series representations, 
it may be useful also to give here the command for the lowest if-type of the discrete 
series with parameter A = discrete of the group U(p, q). The command is: 

>Inf_lowestKtype(p,q, discrete) 

Example 8 



> Inf _lowestKtype ( [ [5/2 , -3/2] , [3/2 , -5/2] ] , 2 , 2) ; 
[[7/2,-3/2], [3/2, -7/2]] 



Similarly, we may want to parametrize a representation of K by its highest 
weight. Then the command is: 

> voganlowestKtype (discrete ,p,q) 
Example 9 

> voganlowestKtype ( [ [5/2 , -3/2] , [3/2 , -5/2] ] , 2 , 2) ; 
[[3,-1], [1, -3]] 
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Let us finally recall the simple case of the multiplicity function for G = U(2, 1) 
(see Sec. 14,11 for the notations). 

Choose as positive compact root the root e% — e2 and denote by w the corre- 
sponding simple reflection. Fix A = [[Ai, A2], [A3]] the Harish Chandra parameter 
for a discrete series representation. We assume A regular and Ai > A2- There 
are three chambers Ci, C2, C3 and hence three systems of positive roots containing 
ei — e%. Precisely ci corresponds to the positive system {e\ — e-ii^-X — — 63}, 
C2 corresponds to the positive system {e\ — e2,ei — e3,e3 — £2) and C3 to {e± — 
e2,e 3 - ei,e 3 - e 2 }. 

We examine the situation in which the discrete series parameter belongs to one 
of these chambers. FigUJand Figj2j picture the two chambers ci, C2 and evidentiate 
the values for m^, when A is in the chamber. The black lines mark the chambers 
containing the compact root e\ — e<i and the red ones the system of positive roots 
for the given chamber. 




Figure 1: m x for the chamber Ci of U(2,l). 
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Figure 2: for the chamber C2 of U(2,l) 



We give some examples concerning the two situations. The parameter A in 
FigUJis 

hol21 := [[2,1], [-3]] 
In Figure [2] the parameter A is 

aba := [[2,-3], [1]]. 
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2 Mathematical background 



2.1 Notations 



G, q semi-simple real Lie group, Lie algebra of G. 

K, t maximal compact subgroup of G, Lie algebra of K. 

T, t maximal torus of K, Lie algebra of T. 

Pet* lattice of weights of T. 

A+, A + (A) C t* system of positive roots for G. 

A+, A+(A) system of positive compact roots. 

A+, A+(A) system of positive non compact roots. 

A + (^4, B) system of positive roots of parabolic type determined by (A, B). 

a, a c positive chamber for G and K respectively 

p, pf,Pn half the sum of positive, positive compact, positive non compact roots. 

P ,P{ set of G admissible, K admissible parameters. 

Pg,P£ set of G admissible and regular , K admissible and regular parameters. 

U r-dimensional real vector space; x € U. 

V dual vector space of U, h G V. 

( , ) the pairing between U and V. 

Vz lattice of V. 

Uz dual lattice in U. 

A + a sequence of vectors in V%; a € A + . 

t a tope. 

T torus U/Ur, t € T. 

F finite subset of T. 

II_4+ (h) polytope defined by A + . 

Na+ (h) number of integral points for H"_4+ (h). 

c chamber. 

JK C Jeffrey-Kirwan residue. 

Ires iterated residue. 



2.2 Blattner formula and multiplicities 

Let G be a reductive connected linear Lie group with Lie algebra and denote by 
K a maximal compact subgroup of G with Lie algebra t. 

We assume that the ranks of G and K are equal. Under this hypothesis the 
group G has discrete series representations. Recall Harish-Chandra's parametriza- 
tion of discrete series representations. We choose a compact Cartan subgroup 
T C K with Lie algebra t. Let P C t* be the lattice of weights of T. They 
correspond to characters of T. Here if A € P, the corresponding character of T 
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is e lX . Let A + C P be a positive system of roots and P = \ Sa!eA+ a - Then the 
subset p + P C t* does not depend of the choice of the positive system A + . We 
denote it by P s . We denote by P£ C P 8 the subset of ej-regular elements. We 
can similarly define P { r . We denote by W c the Weyl group of K. For any A E Pg, 
Harish-Chandra defined a discrete series representation 7r\ Elements of P£ are 
called Harish-Chandra parameters for G. Two representations, 7r A and ir x , co- 
incide precisely when their parameters A, A' are related by an an element of W c . 
Thus the set of discrete series representations is parametrized by P^/W c . 

In the same way we can parametrize the set K of classes of irreducible finite 
dimensional representations of K by their Harish-Chandra parameter p £ P£ /W c - 
Once a positive system of compact roots is chosen, an element p € P£ can be 
conjugated to a unique regular element in the corresponding positive chamber 
o c C t* for the compact roots. We denote by t u € K, or simply by p E K, the 
corresponding representation. 

A discrete series representation ir x is i^-finite: 

7rA Ik = Yl m W- 

To determine m x , that is the multiplicity of the K-type r M in 7r A , is a basic problem 
in representation theory. 

Blattner's formula, |14j . gives an answer to this problem. We need to introduce 
a little more notations before stating it. 

We let A be the root system for q with respect to t, A c the system of compact 
roots , that is the roots of t with respect to t, and A n the system of noncompact 
roots. We let A + be the unique positive system for A with respect to which A 
is dominant. We write p = \ J2atA+ a > Pc = \ S ae A+ a and P™ = \ ^aeA+ a 
where A+ = A + n A c and A+ = A+ n A n . Therefore P g = p + P, P t = p c + P 
and if £ € P s , then £ + p n € P(. 

Write o, a c C t* for the (closed) positive chambers corresponding to A + and 
A+. We will also write A+(A), A+(A) and A+(A) for A+, A+, A+ if necessary to 
stress that these systems depend on A. 
Then for p € P^ Pi o c , Blattner's formula says: 

(1) m X = ^2 e { w )^n{wp - A - p n ) 

where, given 7 G t*, we define V n {^) to be the number of distinct ways in which 7 
can be written as a sum of positive noncompact roots (recall our identification 
for which A, A c C t*). The number V n {l) is a well-defined integer, since the 
elements of A+ span a cone which contains no straight lines. As usual, e(w) will 
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stand for the sign of w. Remark that, as p + p c and A + p n + p c are weights of T, 
the element wp — A — p n is a weight of T. 

It is convenient to extend the definition of to an antisymmetric function on 
PjT. As we observed already an element p £ P£ can be conjugated to a unique 
regular element in the corresponding positive chamber a c C t* for compact roots, 
via an element w € W c . Thus we define = e(w)m^ lfl . Of course with this 
generalization the multiplicity of the K-type p is |m^| and we can complete our 
picture for U(2, 1), Figj3l in the following way: 




Figure 3: as antisymmetric function on U(2, 1). 

The representation Ti owest with Harish-Chandra parameter pi owes t = A + p n is the 
lowest if- type of the representation ir x and occurs with multiplicity 1. It is, in 
general, difficult to compute for general p. 

We will use Blattner's formula to compute m^. Our algorithm is based on a 
general scheme for computing partition functions using multidimensional residues. 
Note that the presence of signs in Blattner formula doesn't even allow to say if a 
K -type appears without fully computing its multiplicity. 

Recall that if A + is a positive root system of a semi-simple Lie algebra, the 
formula for the partition function has been used to compute tensor product de- 
composition or weight multiplicities ([1],[6]). 
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2.3 Polynomial behavior of the Duistermaat-Heckman 
measure 

Recall Paradan's results ([E]) on the behavior of the function m x and its support. 
For this it is useful to first recall the semi-classical analog of Blattner formula. 

Let us fix a positive system of roots A + and consider the corresponding (closed) 
positive chambers a, a c C t* for A + and A+. Our parameter A varies in a and it is 
non singular. In this subsection, the integrality condition A G PJ* is not required. 

Let 0\ C g* be the coadjoint orbit of A. It is a symplectic manifold, and thus 
is provided with a Liouville measure d/3\. Let p : 0\ — > t* be the projection. This 
is a proper map. Each coadjoint i^-orbit in t* intersect a c C t* C f . Thus the 
projection of 0\ on t* is entirely determined by its intersection with a c . We recall 
that the set p(0\) n a c is a closed convex polyhedron. 

Definition 10 The Kirwan polyhedron Kirwan(A) is the polyhedron p(0\) D a c . 

As far as we know, there are no algorithm to determine the Kirwan polyhedron. 

A weak result on the support of Kirwan(A) is that Kirwan(A) is contained in 
A + Cone(A+) where Cone(A+) is the cone generated by positive non compact 
roots. 

The push- forward of the measure d/3\ along the projection p : 0\ — > t* gives 
us an invariant positive measure on t*. By quotienting this measure by the signed 
Liouville measures df^Kp, of the coadjoint orbits K[i in t*, we obtain a W c -anti- 
invariant measure dF x on t* . More precisely, for <fi a test function on t* , 

(2) j Q mWW) = j£ dF\^ dft^WC/)) • 

Here is the locally constant function on t* anti-invariant by W c and equal to 1 
on the interior of the positive chamber a c . We refer to dF x as the Duistermaat- 
Heckman measure. 

If A + = [ai, . . . , a^] is a sequence of elements in t* spanning a pointed cone, 
the multispline distribution is defined by the following formula. For <f) a test 
function on t*: 

roo poo N 

(3) (Y A+ ,<f>)= ••• / ^{S^Ua^dtf-dtN. 

Jo Jo i=1 

Then, for A G o and p, G t* , we have the following result due to Duflo-Heckman- 
Vergne. [15] : 

(4) dF x {^)= Y, <w)w(6 x *Y At ). 
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where A+ is the system of positive noncompact roots defined by o. 

To simplify our next statements, assume that G is semi-simple and has no 
compact factors. Then t* is generated by non compact roots. Recall that the 
spline function is given by a locally polynomial function on t* well defined outside 
a finite number of hyperplanes (see the description later). Choosing the Lebesgue 
measure dh associated to the root lattice, we may identify the measure Y A + to a 
function, denoted by Y^- Similarly, we identify the measure dF x to a function F x 
on t*. Then we have, almost everywhere, the semi-classical analogue of Blattner 
formula: 

F x (li) = e(w)Y+( Wt x-X) 

wew c 

A E a and regular. 

The (anti- invariant) function F x (li) restricted to the positive compact chamber 
o c is a non negative measure with support the Kirwan polyhedron Kirwan(A). 

It follows from the study of spline functions that there exists a finite number of 
open polyhedral cones R 1 in a x a c (so that the union of the cones R^ cover a x a c ) 
and polynomial functions p 1 on a x a c such that F x ([i) is given, for A E o, ll E 
o c , (A, li) E R\ by the polynomial p l (X,ii) on R l (X) = {ll E o c , (A,//) E R 1 }. In 
particular the Kirwan polyhedron Kirwan(A) is the union of the regions i? l (A) for 
which the polynomial p % restricted to R l (X) is not equal to 0. In fact the functions 
p l (X,Li) are linear combinations of polynomial functions of uiX — ll where w are 
some elements of W c . 

If R is an open cone in o x a c such that F x (ll) is given by a polynomial formula 
p R (X,/j.) when (A,/i) G R, X € a, we say that R is a domain of polynomiality and 
that p R is the local polynomial for F x on R. 

Let us finally recall that the local polynomials p R belong to some particular 
space of polynomials satisfying some system of partial differential equations. For a 
a non compact root, consider the derivative d a . We say that an hyperplane H € t* 
is admissible for A n if H is spanned by a subset of dimt — 1 non compact roots, 
that is roots in A n . We denote by H n the set of admissible hyperplanes for A n . 

Definition 11 A polynomial p on t* is in the Dahmen-Micchelli space D(A+) if 
p satisfies the system of equations: 

( n 9 <*)p= q 

aeA+\Q 

for any Q G U n - 
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Remark that the space Z)(A+) depends only of A n and not of a choice of A+. 
Then, the following result follows from Dahmen-Micchelli theory of the splines. 

Proposition 12 For any domain of polynomiality R, the polynomial fi — > p R (X, /i) 
belongs to the space D(A+). 

2.4 Quasi-polynomiality results 

Let us come back to the discrete setting. Let us fix as before a positive system of 
roots A + and consider the corresponding chambers a, a c C t* for A + and A+. Fix 
A € Pg n a and fi € P£ H a c . We can then define m^, the multiplicity of r M in the 
discrete series ir x . 

By definition, a quasipolynomial function on a lattice I is a function on L 
which coincides with a polynomial on each coset of some sublattice L' of finite 
index in L. The subsets P s ,Pt are shifted lattices and we may say that a function 
k on P t r is quasi polynomial on P g x Pi if the shifted function k(X — p, fi — p c ) is 
quasipolynomial on the lattice P x P. 

Theorem 13 Let R be a domain of polynomiality in a x o c for the Duistermaat- 
Heckman measure. Then there exists a quasi polynomial function P R on P g x P^ 
such that = P R (X, n) for any (A, /i) 6 i?fl (P g x P[), A G a, [i G a c . 

(In fact the functions P R are linear combinations of quasi polynomial functions 
of w\ — \x where w are some elements of W c .) 

The iT-types occurring with non zero multiplicity in ir x are such that /u is in 
the interior of the Kirwan polyhedron Kirwan(A). In particular the lowest K-type 
Slowest is i n the interior of Kirwan(A). In particular all the .ff -types occurring 
with non zero multiplicity in 7r A are such that fi is in the interior of the cone 
A + Cone(A+). We believe they are contained in the cone Slowest + Cone(A+), 
but we do not know if this assertion is true or not (by Vogan's theorem, they are 
contained in Slowest + Cone(A + )) . 

If v & i*, we say that v is an asymptotic direction, if the line Slowest + tv 
is contained in Kirwan(A) for all t > 0. The set of asymptotic directions form a 
cone, which determines the wave-front set of tt )k \k- 

For the holomorphic discrete series, the descrition of the cone of asymptotic 
diections is known. In fact if the lowest i<C-type of 7r A is a one dimensional rep- 
resentation of K, the exact support of the function m x has been determined by 
Schmid. 

We will explain in Section 13.71 how to compute regions of polynomiality R and 
the quasi-polynomial P R . 
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The quasi polynomials P R are in some particular space of quasi polynomials 
satisfying some system of partial difference equations. For a a non compact root, 
consider the difference operator V a acting on Z valued functions on Pf by 

(V a k)(n) = k(fi) - k(n - a). 

Definition 14 A quasi polynomial L on P% is in the Dahmen-Micchelli space 
DM(A+) if p satisfies the system of equations: 

( n v q )l=o 

for any Q € U n . 

Then, the following result follows from Dahmen-Micchelli theory of partition func- 
tions. 

Proposition 15 The quasi polynomial \i — > P R (X, /x) belongs to the space DM(AJ). 

2.5 Aim of the algorithm: what can we do? 

Our algorithm addresses the following questions for U(p, q). All of these questions 
will be analyzed in more details in Sec J3.71 

2.5.1 Numeric 

We enter as input two parameters A, /i € PJ* x P { r . The output is the integer m^, 
see Sec l3.7,ll 

2.5.2 Regions of polynomiality 

The input is two parameters Ao, (J>o € P r x P£. Let o, a c be the chambers determined 
by Ao and fj,Q. We also give two symbolic parameters A, jj,. 

Then the output is a closed cone R(\q,iiq) C o © a c described by linear in- 
equations in A, fj,, containing (Ao,//o) an d a quasi-polynomial P in (A,^) such that 
ml = P(A, n) for any (A, /i) G R(X , Mo ) n (P; x P{). 

We worked out part of this program for U(p,q), but it is still not fully imple- 
mented. 

In particular, for the moment, we are not able to produce a cover of a x a c 
by such regions. The number of regions needed grows very fast with the rank. 
Furthermore, we are not able to decide when we have finished to cover a x o c . 
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2.5.3 Asymptotic directions 

We implemented (for U(p, q)) a simpler question which gives a test for asymptotic 
directions. 

Let's consider as input parameters Ao in P g r and a weight v G a c . Let //q be the 
lowest if -type of 7r Ao . The line t i— > (Ao, /^o + ^) cross domains of polynomiality 
R 1 at a certain finite number of points < t\ < t% < ■ ■ ■ < t s . Let us define 
to = 0,i s +i = oo. Then we study the function P(t) = m^° +to , t G N,t > 0. 

We can find polynomials q\t u t i+1 ] of degree bounded by pq — (p + q + 1) such 
that P(i) = tj,]! (t) when t{ < t < ti+i for i = 0, . . . , s and f G N. In particular, 
the direction v belongs to the asymptotic cone of the Kirwan polyhedron, if and 
only if our last polynomial q\t a CO ] is non zero, see Sec J3.7.3l 

As we discussed in the first part, if the intervals are two small, these polyno- 
mials are not uniquely determined. However the last interval is infinite, and the 
last polynomial is well determined. If this last data is non zero, then v is in the 
wave front set of tt x . The reciproc is not entirely clear. Indeed for a direction to 
be in the wave front set, it is sufficient to be approached in the projective space 
by lines M + /i n with \i n G P£ PI a c such that the multiplicity mj n is non zero, and 
the sequence /i n is going to the infinity in a c . Thus we do not know if a rational 
line fiQ + tv contained in the Kirwan polytope could totally avoid the support of 
the function We do not think this is possible. 

3 Partition functions: the general scheme. 
3.1 Definitions 

Let U be a r-dimensional real vector space and V be its dual vector space. We 
fix the choice of a Lebesgue measure dh on V. Consider a list A + of non-zero 
generators for V given by 

A + = [ai,a 2 , ■ ■ ■ ,a N }. 

We recall several results concerning partition functions that appear in [3] in 
the general context. However, let us describe right away the system of vectors 
A + (A,B) C A r that will appear in our programs and that describe parabolic 
subsystem of A r , (as we said the same method could be applied to other parabolic 
root systems). 

Example 16 • Let E be an r-dimensional vector space with basis e« (i = 1, 
. . . , r + 1). Consider the sequence 

A+ = [e i -e j \l<i<j <r+ 1]. 
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This is a system of positive roots of type A r . We let V to be the vector space 

C r+1 r+1 ^| 

v = h=j2 h * e * e E I Yl hi = [ • 

l i=l i=l ) 

Let V% be the lattice spanned by . We may identify V with W by h ^ 
\h\ , /t2 , . . ■ , h r ] . In this identification the lattice V% is identified with HI . 

• Let A, B be two complementary subsets of [1, 2, 3, . . . ,p + q] with \A\ = p and 
\B\ = q. Let r = p + q — 1. We define A + (^4, B) as the sublist of defined 
by 

A + (A, B) = [a - ej 1 1 < i < j < r + 1], with i G A,j G B or i G B,j G A. 

These systems are the system of positive roots for the maximal parabolics of 
g[(r + 1), for different choices of orders. 

Let us go back to the general scheme. 
For any subset S of V, we denote by C(S) the convex cone generated by non- 
negative linear combinations of elements of S. We assume that the convex cone 
C(A + ) is acute in V with non-empty interior. 

If S is a subset of V, we denote by < S > the vector space spanned by S. 

Definition 17 A hyperplane H inV is ^4 + -admissible if it is spanned by a set of 
vectors of A + . 

When A + = A^(A) or A + (A, B) then an A + -admissible hyperplane will be 
also called a noncompact wall. 

Chambers 

Let V s i n g(A + ) be the union of the boundaries of the cones C(S), where S ranges 
over all the subsets of A + . The complement of V s i ng (A + ) in V is by definition the 
open set C reg (A + ) of regular elements. A connected component c of C reg {A + ) is 
called a chamber of C{A + ). Remark that, in our definition, the complement of the 
cone C(A + ) in V is a chamber that we call the exterior chamber. The chambers 
contained in C(A + ), that we will call interior chambers, are open convex cones. 
Sometimes chambers are called cells or big cells by other authors. 
The faces of the interior chambers span admissible hyperplanes. 
The following pictures illustrate the situation for the (interior) chambers in the 
case of A^ , FigUl and the (interior) chambers for various subsystems of A^ , of 
type A+(A,B) relative to 17(2,2), FigEJ 
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e 2 - e 3 




e, - e 2 e 3 - e 4 

A + ([1 , 3], [2, 4]) 



Figure 5: Parabolic subsystems of A% 
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Polytopes 

We consider the space HH N with its standard basis cjj and Lebesgue measure dx. 
If x = YliLi x i UJ i R N i we simply write x = (x±, . . . ,Xjy). Consider the surjective 
map A : R N -> V defined by Afa) = a*. 

If h G V, we define the convex polytope II4+ (/t) consisting of all non-negative 
solutions of the system of r linear equations Yli=i x i a i = h that is 

U A +(h) = {x = (xi,...,x N ) eR N \Ax = h, x t > 0} . 

We call H A +(h) a partition polytope (associated to A + and /i). 

We identify the spline distribution (by Formula [3|) to a function still 
denoted by Y4+ using dh. 

Recall the following theorem, which follows right away from Fubini theorem, 

(CD!, HBO 

Theorem 18 The value of the spline function Y4+ at h is the volume of the 
partition polytope U^+(h) for the quotient measure dx/dh. 

The spline function Y4+ is given by a polynomial formula on each interior 
chamber. It is identically equal to on the exterior chamber. 

Partition functions 

Let Vz be a lattice in V and suppose now that the elements oti of our sequence 
A + belong to the lattice V%. If h G Vz we define N A +(h) = \U A +{h) C\Z N \, the 
number of integral points in the partition polytope U A +(h). 

Thus N A +(h) is the number of solutions (x±,X2, ...,xn), in non-negative in- 
tegers Xj, of the equation Ylj=i x j a j = h- 

The function h *— > N A +{h) is called the partition function of A + . We refer to it as 
Kostant partition function. 

We will see after stating Theorem [23] that h h-> N a + (h) is quasipolynomial on 
each chamber. 

Let us recall briefly the theory that allows to compute Kostant partition func- 
tions. 

JefFrey-Kirwan residue 

Let v be a subset of {1, 2, . . . , N}. We will say that v is generating (respectively 
basic) if the set {aj | i € u} generates (respectively is a basis of) the vector space 
V. We write Bases(*4 + ) for the set of basic subsets. 

Let 1Z A + be the ring of rational functions on U, the dual vector space to V, 
with poles on hyperplanes determined by kernel of elements a 6 A + . 
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7£_4+ is Z-graded by degree. Every function in TZj^+ of degree — r decomposes 
(see [5]) as the sum of basic fractions f a , f a = n 1 , a € Bases(*4 + ) and 
degenerate fractions; here degenerate fractions are those for which the linear forms 
in the denominator do not span V. 

Now having fixed a chamber c, we define a functional ,lK c (f a ) on TZ^+ called 
the Jeffrey-Kirwan residue (or JK residue) as follows: 



(5) JK c (/ ff ) 




- 1 , ifcCC(a), 
if cn C(<t) = 



where a 6 Bases(^4 + ) and vol(cr) is the volume of the parallelotope ^ =1 [0, l]aj 
computed for the measure dh. 

There exists a linear form JK C , that we call the Jeffrey-Kirwan residue, on 1Z^+ 
such that JK C takes the above values on the elements f a , and is equal to on a 
degenerate fraction or on a rational function of pure degree different from — r. 

If c is the exterior chamber, then clearly JK C is equal to 0, as c is not contained 
in C(A+). 

We may go further and extend the definition of the Jeffrey-Kirwan residue to 
the space TZ^ which is the space consisting of functions Pj Q where Q is a product 
of powers of the linear forms oti and P = Yl^o is a formal power series. Then 
we just define, if Q is of degree q, 

JK C (P/Q) = JK c (P 9 _ r /Q) 

as the JK residue of the component of degree — r of P/Q. 

3.2 Spline functions and Kostant partition function 

Let us recall the formulae for the spline function and for Nj L +{h). 

Definition 19 Let c be an chamber contained in the cone C{A + ). Define the 
function Y c on V by 



Y c A+ {h) = JK C 



e h 



niV I 

More explicitly, as JK C vanishes outside the degree — r, we have 

1 / h N ~ r \ 



We thus see (h) is an homogeneous polynomial on V. 
The proof of the following theorem is immediate (|10|). 
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Theorem 20 Let Y A +(h) be the multispline function associated to A + . Let c be 
a chamber contained in the cone C(A + ). 
We have for h G c 

Y A+ (h)=Y A+ (h). 

Remark 21 According to Theorem \18\, this theorem gives the formula for the 
volume V A +(h) of the partition polytope U A +(h). 

Let us now give the residue formula for the number of integral points N A + of 
the partition polytope H A +(h). 

Consider the torus T = U/U% where U is the dual vector space to V and 
Uz C U is the dual lattice to V%. If G G U, we denote by g its image in T. 

For a G Bases(^4 + ) we denote by T(a) the subset of T defined by 



T(o) = {geT efoW^TG) = i for alia € a, G a representative of g G U/Uz\ . 
The set T(a) is a finite subset of T. 

For G G U and h G V, consider the Kostant function K(G, h) on U defined 

by 

Ah,2/Ky/-lG+u) 

(6) K(G, h)(u) = 



Remark 22 If h G Vz, the function K(G,h) depends only of the class g of G in 
U/U z . 

The function K(G, h)(u) is an element of TZ A . 

Indeed if we write 1(g) = ji 1 < i < N, e ~^^^G) _ A ; then 

( 7 ) k(g, „)(«) = ^^» ^m. 

where ip 9 (u) is the holomorphic function of u (in a neighborhood of zero) defined 

by 

ih9( u ) = TT x TT - 

By taking the Taylor series of e( h > u 'if) 9 (u) at u = 0, we see that the function 
ii — > K(G, h)(u) on ?7 defines an element of 1Z A . If c is a chamber of C(A + ), the 
Jeffrey-Kirwan residue JK c (K(g, h)) is thus well defined. 
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Definition 23 Let c be a chamber. Let F be a finite subset of U . We define the 
function on V by 

N^O) = vol(V/V z ,dh) £ JK t (K(G,h)) 

G&F 

where vol(V/Vz, dh) is the volume of the fundamental domain ofVz for dh. 

Finally introduce the zonotope Z(A + ) to be the convex polyhedra defined by 

N 
i=l 

When A + is fixed, we just write Z = Z(A + ), and if C is a set, we denote by C — Z 
the set of elements {£ — z} where £ € C and z € Z. 

The following theorem is due to Szenes-Vergne [3]. It generalizes [Tj, [13] and 

0- 

Theorem 24 Let c be a chamber. Let F be a finite subset of U . Assume that for 
any a £ Bases(^4 + ) such that c C C(a), we have T{a) C F/Uz- 
Then for h G V% H (c — Z), we have 

N A+ (h)=N%(h). 

We choose for any chamber c such a finite set F such that all elements g € F/Uz 
have finite order and such that F satisfies the condition: 

(C) for any a € Bases(^l + ) such that c C C(a), we have T(<r) C F/Uz- 

It is possible to achieve this, for example choosing a set F of representatives 
of -U% modulo Uz, where p is that pUz is contained in J2iea f° r an y basis a. 

We now simply denote N C ' F by N c , leaving implicit the choice of the finite set 

F. 

Remark 25 • Observe that N^ + (/i) does not depend on the measure dh, as it 
should be. 

• If c is the exterior chamber, then N^ + (/i) = 0. In our algorithm, we are not 
knowing in advance if the point h belongs to the cone C(A + ) or not, so that this 
remark is not as stupid as it looks. 

• Observe also that if c is an interior chamber, then c — Z contains the closure c 
of c, while if c is the exterior chamber c — Z = c. For an interior chamber, usually 
the set c — Z intersected with the lattice Vz is strictly larger than c intersected with 
V%. This fact will be important for computing shifted partition functions, as we 
will explain later. 
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Let us explain the behavior of the partition function N A + on the domain c — Z. 
We first explain the case of an unimodular system. 

Definition 26 The system A + is unimodular if each a £ Bases(.4 + ) is a Z-basis 
ofV z . 

Example 27 It is easy to see that is unimodular, so is any subsystem. 

Thus if A + is unimodular, the set F = {0} satisfies the condition (C) and we 
choose this set F. 

Proposition 28 If A + is unimodular, the function (h) is a polynomial func- 
tion on V. 



Proof. We have just to consider K(G, h) = K(0, h) and we can write 



K(0,h)(u) 



,(h,u) Jh,u) nJV 



nf=i(i - e-<*.«>) nf=i (<*,«> nf=i(i - e-<*.«>) 



where Fpftp = 12k=o ^k( u ) is a holomorphic function of u in a neighbor- 

hood of with ipo(u) = 1. 

It follows that N^ + (/i) is given by the following polynomial function of h 

( e^ \ 
N A+ (h) = vo\(V/Vz,dh)JK c i-^- -x 

VIli=i <«<.«> fc =o / 

Note that the function is a polynomial function of degree N — r whose 

homogeneous component of degree N—r is the function (h), that is the volume 
of the polytope. 

Let us now consider the general case where F is no longer reduced to {0}. 
For example for parabolic root systems of B r , C r , D r , the set F satisfying the 
condition (C) cannot longer be taken as equal to {0}. 

We recall that an exponential polynomial function is a linear combination of 
exponential functions multiplied by polynomials. 

Proposition 29 The function (h) is an exponential polynomial function on 
V and the restriction o/N^ + (/i) to Vz is a quasipolynomial function on Vz- 
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Proof. Let us denote by ip 9 (u) = Y0k=o the ser i es development of the 

holomorphic function ip 9 appearing in formula ([7]). Then we see that JK C (K(G, h)) 
equals 



e {h,u) 



(9) e^^^JK -r{u) 

\ \UaI{g){^ U ) J 

>— s |/( ^l~ r 1 ( (h v \\^3)\-r-k 



The function 



» JK C A^L^ ^ (u) 



is a polynomial function of /i of degree |X(c/) | — t — Ar. Thus we see that ,1K C (K(G, h)) 
is the product of the exponential function e ( h ' 27T V-iG) a p iy nom ial function of 
h. _ 

Furthermore, if g is of order p and h varies in Vz, the function h *— > e (' l . 27r v /3 TG) 
is constant on each coset h + pVz of the lattice pVz- 

Return to the computation of the partition function N^+(h). Thus we see 
that when h varies in (c — Z) n Vz, we have that N^+(h) coincide with the quasi 
polynomial function N^ + (/i) above. Note that its highest degree component is 
polynomial and is again the function (h), the volume of the polytope II4+ (h). 

The quasipolynomial nature of the integral-point counting functions iV^ stems 
precisely from the root of unity in formula ([9]). 

Furthermore for parabolic root systems of type B, C, and D, these roots of 
unity are of order 2, as in the following example. Thus we summarize the properties 
of our partition functions in the following remark: 

Remark 30 • is unimodular, that is we can choose F = in Theorem 
\24\ and thus the partition function N$ for any subset <3> of coincide with 
a polynomial function on each domain c — Z(Q). 

• The integral-point counting functions iV$ for any subsystem of B r ,C r ,D r 
coincide with quasipolynomials with period 2 on each domain c — Z(Q). 

We now compute the number of integral points in two different situations: a 
non unimodular case and a unimodular one. We treat the non unimodular case 
first. 

Example 31 Here V is a vector space with real coordinates and basis e\,C2 and 
U = V* has dual basis e 1 , e 2 . We write v = Y^a=i v i e i e ^ an ^ u = Yli=i ^* e * ^ ^ 
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for elements in V and U respectively. Let us compute the number of integral points 
for the positive non compact root system occuring for the holomorphic discrete 
series of SO(5, C) : that is we fix A + := {e±, e 2 , &i + &2, &i — £2} and A + = A+ := 
{ei,ei + e2,ei — e 2 }. We also write a vector h = h\e\ + h 2 e 2 in the cone C(A + ) 
as (hi,h 2 ). Of course, the calculation can be done by hand, but we illustrate the 
method in this very simple example. 

Observe that the root lattice is 7Le\®7Le 2 andvol (V/Vz, dh) = 1 for the measure 
dh = dh\dh 2 . 

There are two chambers, namely t\ = C{{e\ + e 2 , ei}) and C2 = C({ei, e\ — 62}). 
Now let us compute the Jeffrey-Kirwan residues on the chambers. 
We have for example: 



JK 



JK C1 ( , 1 . A = 1, JK C „ ( , \ v ) = 0, 

c i \ui(ui+u 2 ) I C2 \Ul(u 1 +U2) I 

( 1 \ _ 1 tk ( 1 \ - I 

c i \{u 1 +u 2 ){u 1 -u 2 ) I 2> JJ:Vc 2 \{u 1 +u 2 ){u 1 -u 2 ) I 2' 

JK C1 ( —r± A = 0, JK C9 ( —r^ A = 1 



For i/ie number of integral points, we first note that F = {(0,0), (1/2, 1/2)}. 
Consequently N A +(h) is equal to the Jeffrey-Kirwan residue of fi = K((l,l), h) 

plus f 2 = if ((1/2, 1/2), h). We rewrite the series fj (j = 1, 2) as fj = /■ x 

e Ulhl+U2h2 /mim + u 2 )(u 1 -u 2 ) where 



, Ml Ml + M2 Ml - M 2 

/l 



1 — e~ ui 1 — e _ ("i+"2) 1 — e -("i-«2) ' 

/ 2 = — M!+M2 M!-M 2 { _ 1)hl+ H 2 

Using the series expansions l J L e - x = l+^x+ ^x 2 +0(x 3 ) and = ^x+0(x 2 ), 

we obtain that the number of integral points is the JK residue of 

fti + 3 , I \ h 2 u 2 



(Ml — M 2 )(m1 + M2) Ml(«l+M 2 ) Ml(Ml — M 2 ) Ml (mi — M 2 ) (mi + M 2 ) (m X + M 2 ) (mi — M 2 ) 
_ fel + | + |(-l)^+^ | § | § h 2 | h 2 

(mi — m 2 )(mi + m 2 ) mi(mi+m 2 ) mi(mi — m 2 ) mi(mi+m 2 ) (mi + m 2 )(mi — m 2 ) 
We then obtain: 



N A +(h) = ih 1 + i(-i)''i+h2 + |_i h2) i/ft 6Cl , 



iVoie that the functions N A + agree on walls, that is h 2 = 0, and the formulae 
above are valid on the closures of the chambers. 
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The second example treats the unimodular case of A+, see Example 1161 Since 
we have identified V with W , then we have a canonical identification of U = V* 
with R r defined by duality: u € W to u = Yll=i u i e% ^ where e 1 is the dual 
basis to ej. Thus the root ej — ej (1 < i < j < r) produces the linear function 
m — uj on U, while the root ej — e r+ i produces the linear function Uj. Recall also 
the identification h = ^i=i ^« e * = [hi, ■ ■ ■ ,h r ], 

We compute the number of integral points for the parabolic subsystems of 
U(2, 2) illustrated in FigEJ 

Example 32 We consider the 3 different systems of non compact roots as de- 
scribed in Fig\f^ and give the formulae for the partition function. 

1. J/A+ = A+([1,4],[2,3]) then 

h± + h 2 + 1 if h £ Ci, 

hi + h-2 + h-3 + 1 if h£c2, 

hi + hz + 1 if h e C3, 

hi + 1 if h £ C4 




5. I/A+ = A+([l,2],[3,4]) to 



iV A+ (ft) 




«/ /i 6 ci, 

if h£ c 2 , 

if he C3, 

i/ 7i 6 C4 



5. //A+ = A+([1,3],[2,4]) to 

r 1 + fti + A 2 if he a 

N A+ (h) = i l + hi+h,2 + h 3 if h e c 2 , 
[ 1 + if hec 3 , 

We have to compute the Jeffrey-Kirwan residue of the function 

computation is immediate since we need only term of degree one for the expansion 
°f fi ■ We omit the details. Remark though that once again the formulae agree on 
walls as it should be. 

3.3 Shifted partition functions 

Let us consider as before our lattice V% and our sequence A + of elements of Vz- 
Let 

Pn = 2 ^ a ' 

aeA+ 
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We introduce 



P n = Pn + V z . 



Thus for any p, £ P n , the function (p — p n ) is well defined. 

Let H be the complement of all admissible hyperplanes, that is hyperplanes 
generated by elements of A + , DeflTTl 

Definition 33 A tope is a connected component of the open subset V — H ofV. 

We choose once for all a finite set F of elements G of U, so that the image of 
elements g cover all groups T{a). 

If t is a tope, then r is contained in a unique chamber c, and we denote by 
the exponential polynomial function given in Definition 1231 If r is not 
contained in C(A + ), then = 0. 

The closures of the topes r form a cover of V. A consequence of Theorem 1241 
is the following. 

Theorem 34 For any tope r such that /i € r n P n , we have 



3.4 A formula for the Jeffrey-Kirwan residue 

Having stated a formula for partition functions (or shifted partition functions) 
in terms of JK C , we will explicit it using the notion of maximal proper nested 
sets, as developed in [12], and the notion of iterated residues. The algorithmic 
implementation of this formula is working in a quite impressive way, at least for 
low dimension. 

This general scheme will be then be applied to Blattners' formula. 
3.4.1 Iterated residue 

If / is a meromorphic function of one variable z with a pole of order less than or 
equal to k at z = 0, then we can write f(z) = Q(z)/z k , where Q(z) is a holomorphic 
function near z = 0. If the Taylor series of Q is given by Q(z) = Yl^LoQsZ s , then 
as usual the residue at z = of the function f(z) = QsZ s ~ h is the coefficient 

of 1/z, that is, qk-i- We will denote it by res z= o/(z). To compute this residue we 
can either expand Q into a power series and search for the coefficient of z^ 1 , or 
employ the formula 



N A+ (A 4 ~ Pn) = (p - p n ). 



(10) 



res z=0 f(z) 



1 



(d^- 1 (z k f(z)) 



(fc-1)! 



2 = 
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We now introduce the notion of iterated residue on the space TZ^+ . 

Let v = [ai, Q2, • • • , a r ] be an ordered basis of V consisting of elements of 
A + (here we have implicitly renumbered the elements of A + in order that the 
elements of our basis are listed first). We choose a system of coordinates on U 
such that ai(u) = Ui. A function (j) £ TZ^+ is thus written as a rational fraction 
(f)(ui,U2, . . . ,u r ) = Q("|'^'"''" 7 'j where the denominator Q is a product of linear 
forms. 

Definition 35 If <f) £ TZa> the iterated residue IreSj?((/>) of 4> for V is the scalar 

Ires^O) = res Ur . =0 res Ur _ 1=0 • • • res Ul=0 <f)(u 1 ,U2, ...,u r ) 

where each residue is taken assuming that the variables with higher indices are 
considered constants. 

Keep in mind that at each step the residue operation augments the homoge- 
neous degree of a rational function by +1 (as for example res x= o(l / xy) = 1/y) 
so that the iterated residue vanishes on homogeneous elements (j) £ 7£_4, if the 
homogeneous degree of 4> is different from —r. 

Observe that the value of Ires,?((/>) depends on the order of z7. For example, for 
/ = l/(x(y - x)) we have res x=0 reSj /=0 (/) = and reSj /=0 res :r= o(/) = 1. 

Remark 36 Choose any basis ji, 72, . . . , 7 r of V such that ® 3 k=1 oij = ©]. =1 7j 
for every 1 < j < r and such that 71 A 72 A • • • A j r = a% A cti A • • • A a r . Then, by 
induction, it is easy to see that for (j) £ T^-a+ 

res ar=0 • • • res ai=o = res 7r=0 • • • res 71=0 </>. 

Thus given an ordered basis, we may modify 02 by 02 + cct\, . . . , with the 
purpose of getting easier computations. 

As for the usual residue, the iterated residue can be expressed as an integral as 
explained in [3]. This fact allows change of variables. 

3.4.2 Maximal proper nested sets adapted to a vector 

We recall briefly the notion of maximal proper nested set, MNPS in short, and 
some of their properties (see p2J). 

A subset S of A + is complete if S = (S) n A + : here recall that (S) is the 
vector space spanned by S. A complete subset S is called reducible if we can find 
a decomposition V = V\ ffi V2 such that S = Si U S2 with Si C Vi and S2 C V%. 
Otherwise S is said to be irreducible. 
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A set M = {I±, I2, . . . , Jfc} of irreducible subsets of A + is called nested if, given 
any subfamily {I\, . . . , I m } of M such that there exists no i, j with ij C Ij, then 
the set I\ U • • • U I m is complete and the elements are the irreducible components 
of I\ U 1-2 U ■ ■ ■ U I m . Then every maximal nested set M, MNS in short, contains 
A + and has exactly r elements. 

We now recall how to construct all maximal nested sets. We may assume that 
A + is irreducible, otherwise just take one of the irreducible components. If M 
is a maximal nested set, the vector space (M \ A + ) is an hyperplane H, thus an 
admissible hyperplane. 

Definition 37 Let H be a A + -admissible hyperplane. A maximal nested set M 
such that (M \ A + ) = H is said attached to H. 

Given M a MNS for A + attached to H, then (M\A + ) is a MNS for Ht~)A + . 
Therefore maximal nested sets for an irreducible set A + can be determined by 
induction over the set of ^4 + -admissible hyperplanes. 

For computing the Jeffrey-Kirwan residue, we only need some particular MNS's. 
Let us briefly review the main ingredients. 

Fix a total order ht on A + . Let M = {Si, S2, ■ ■ ■ , S^} be a set of subsets of 
A + and choose in each Sj the element c<j maximal for the order given by ht. This 
defines a map O from M to A + and we say that M is proper if G(M) = M is a 
basis of V. We denote by V{A + ) the set of MPNS. 

So we have associated to every maximal proper nested set M an ordered basis, 
by sorting the set M = [a%, ai, ■ ■ ■ , a r ] of elements of A + . 

Let v be an element in V not belonging to any admissible hyperplane. 

Definition 38 Define V(v, A + ) to be the set of M G V(A + ) such that v G C{M) = 
C(ai, . . . ,a r ). 

When there is no possibility of confusion we will drop simply write V(v) for 



We are now ready to state the basic formula for our calculations. 

Theorem 39 ( [12J) Let c be a chamber and let v £ c. Then, for (j) £ T^_a+, we 
have 



Let us finally sketch the algorithm to determine V(v,A + ) without going to 
construct all the MNS's. If A + = Af x A\ is reducible, then V(v,A + ) is the 
product of the corresponding sets V(vi,Af). 



V(v,A + ). 




MeV(v,A+) 
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Assume A + is irreducible. Let 9 be the highest root of the system A + (for our 
order ht). We start by constructing all possible .4 + -admissible hyperplanes H for 
which v and 9 are strictly on the same side of H. In particular, the hyperplane H 
does not contain the highest root. 

Then we compute the projected vector projnv on H parallel to 9: v = proj H v+ 
t9, with proj^ v € H and t > and compute A + n H. If Mh is in V{projnv, A + fl 
H), then M = {op(Mh), A + } is in T-^u, A + ). Running through all hyperplanes H, 
for which v and 9 are strictly on the same side of H, we obtain the set 7 3 (f,^4 + ). 
Let us summarize the scheme of the algorithm in Figure [6j Recall that we have 
as input a regular vector v, and as output the list of all MPNS's belonging to 
V(v,A+). 

for each hyperplane H do 

check if v and 8 are on the same side of H 
if not, then skip this hyperplane 
define the projection proj^-u) of v on H along 
write A fl H as the union of its irreducible components I± U • ■ ■ U Ik 
write v as v\ © ■ ■ ■ © according to the previous decomposition 
for each Ij do 

compute all MPNS's for Vj and Ij 

collect all these MPNS's for Vj and Ij 

end of loop running across Ij 's 
collect all MPNS's for the hyperplane H by taking the product of V(Ij, vj) 
end of loop running across H 's 
return the set of all MPNS's for all hyperplanes 



Figure 6: Algorithm for MPNS's computation (general case) 

In our program, we run this algorithm for an element v not in any admissible 
hyperplane, without knowing in advance if v belongs to the cone C{A + ). The 
algorithm returns a non empty set if and only if v belongs to C(A + ). 

3.5 The Kostant function: another formula for subsys- 
tems of 

In this article, we will be using partition functions for lists A + (A, B) described in 
the Example PT6l These lists are sublists of a system of type (with r = p+q — 1). 
In residue calculation, we can use change of variables and thus use a formula for 
which iterated residues will be easier to compute. 

Let us describe this formula. We will describe it for sublists, eventually, with 
multiplicities of a system A+ . We take the notations of Example [TBI 

Let $ be a sequence of vectors generating V and of the form (a — ej), 1 < i < 
j < (t + 1), eventually with multiplicities. Let my (i < j) be the multiplicity of 



33 



the vector e^ — ej in <3? and define tj = nijj+i + • • • + rrij^+i — 1. We recall our 
identification of V with W and of U = V* with lR r defined by duality. In this 
way, as we already observed the root e^ — ej (1 < i < j < r) produces the linear 
function Uj — Uj on U, while the root — e r+ i produces the linear function u; L . 

We are now ready to give another formula for the Kostant function in this 
situation. 

Theorem 40 Let t be a chamber of C{§). Let h = Y^i=i hi&i = [hi,---, h r ], then 
N%{h) = vol (V/Vz, dh) JK C (/* (h) (u) ) , where 



n (i + 



f*(h)(u) 



i=i 



nC" r+1 n (tii-u,-)' 

i=l l<i<jr'<r 



Thus, when /i G Vz n (c - .£($)), we have 
Example 41 • If & = A+ , then 

m=i(i+^) ,ii+r - i 



rii<i<j< r (^ - «j) x n[=i «< ' 

• Let p, q integers such that p + q = r + 1 and <3? 6e £/te system of positive 
noncompact root for defined by $ = {ej — e^, 1 < i < p, p+1 < j < r+1}. 
ITiai zs $ = A+(A,5) wtfi A = [1,... ,p] and B = [p + 1, . . . ,p + q\. 

Then 

U(hi, h 2 ,-.., h r ){u) = 
(1 + nQ^+g- 1 ■■■(! + u p ) h p+"- 1 (l + Mp+i)^ 1 - 1 ■■■(! + u p+q - 1 ) h p+<>- 1 - 1 
(m - Up+i) ■ ■ ■ (ui - • • • (-Up - u p+1 ) • • • (-Up - n p+<? _i)nin 2 ■ ■ ■ u p 

Proof. The function K(0,h)(u) = e^ h ' u ^ / Yl a£ ^(l — e~^ a ' M ^) computed for the 
system is 



e hiui e h 2 u 2 . . . e h r u r 



K(0,h)(u) = - 

Yl (1 _ e -ui}mi,r+i Y\ (1 — e~^ Ui ~ u i">) m <i 
i=l l<i<j<r 

Note that the change of variable 1 + Zi = e Ui preserves the hyperplanes Ui = 
and Ui = Uj and that Z{ = e Ui — 1 leads to dzi = e Ui dui = (1 + Zi)dui- Thus after 
the change of variable we get the required formulae. 
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3.6 Computation of Kostant partition function: gen- 
eral scheme 

3.6.1 Numeric 

We have as input A + a sequence of vectors in our lattice Vz, a vector h G Vz, and 
we want to compute N A +(h). We will compute it by !! 

N A+ (h) = N A+ (h + p n - p n ). 

We mean: Let r be any tope such that h! = h + p n belongs to the closure of r. 
Using Theorem 1341 then 

N A+ (h) = N A+ (ti - p n ) = W A+ (ti - p n ). 

To compute a tope r containing b! , we can move h + p n in any generic direction 

e. 

Here is an outline of the steps needed to compute the number N A + (h) by the 
formula N A +(h) = N T A+ (h). 

Input: a vector h 6 Vz, and A + a sequence of vectors in V%. 
Output: the number N A +(h). 

• Step 1 Compute the Kostant function 

or more generally compute a set F and the functions K(G, h) for G 6 F. 

• Step 2 Find a small vector e so that if h is in Vz, the vector h + p n + e does 
not belong to any admissible hyperplane. Thus the vector h reg = h+p n + e is 
in a unique tope r. The procedure to obtain h reg is called DefVector nc (h). 

• Step 3 Compute the set All := V{h reg ,A + ) as explained in Figj6j 

• Step 4 Compute N T A+ (h) by computing the iterated residues of K(G,h) 

— > 

associated to the various ordered basis M for M varying in the set All. 
That is compute the number 



out := e^' 2 "^ 1 ^ Yl Ires ^( G > h ) 

GeF M&All 



where M is the ordered basis attached to M. 



Then N A + (h) = out 
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3.6.2 Symbolic 

The previous calculation runs with symbolic parameters. If hfix is an element in 
Vz, we might want to find a tope r such that hfix belongs to the closure of r. 
Then 

N A+ (h)=N r A+ (h) 

will be valid whenever h is in the closure of r. Here is the outline of the algorithm. 
Input: hfix is an element in Vz and A + a sequence of vectors. 
Output: A domain D <ZV and an exponential polynomial function P(h) on 

V. 

The domain D is a closed convex cone in V (described by linear inequations) 
such that hfix is in D. The formula P(h) = N A +(h) is valid whenever h € DP\Vz- 

• Step 1 Consider h as a parameter and compute the Kostant function 
K(h){u) as a function of (h,u) given by 

K(h)(u) = K(0,h)(u) = — ^ ' > - — — 

or more generally compute a set F and the functions h){u) for G £ F, 

as function of (h,u). 

• Step 2 Find a small vector e so that if hfix is in Vz, then the vector 
hfix re g := /i/zx + p„ , + e does not belong to any admissible hyperplane. 

Compute the domain D :=T where r is the unique tope r containing hfix reg . 

• Step 3 Compute the set All := V(hfix reg , A + ) as explained in Figj6j 

• Step 4 Compute N^ + (/i) by computing the iterated residues of K(G,h) 

associated to the various ordered basis M for M varying in the set All, here 
h is treated now as a parameter. 

That is we compute 

out := e^^ 1 ^ Ires^(G,/i) 

GeF MeAii 

where M is the ordered basis attached to M. The output out is an expo- 
nential polynomial function P(h) of h and once again we compute 



N A +{h) = P{h) = out, VheDn 



The domain D is a rational polyhedral cone which includes hfix. 

In practice, this works only for small dimensions and when A + is not too big. 
We will program variations of these algorithms, with less ambitious goals. 
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3.7 Computation of Blattner formula: general scheme. 



In this subsection, we summarize the steps to compute Blattner's formula and the 
general scheme to obtain the region of polynomiality. The relative algorithms will 
be outlined in Section [6j 

Let G, K, T be given as in Section 12.21 Let A n C t* be the list of noncompact 
roots. 

Our inputs are A £ P g r and p £ P[. The goal is the study of the function 
p — )• m^. Let A + = A+(A) and recall that in this case a .4 + -admissible hyperplane 
is called a noncompact wall. 

We use Blattner's formula. In our notations: 



3.7.1 Numeric 

Input A £ P r and p £ P[ . 
Output a number. 

The algorithm is clear: 

• Compute A + = A+(A) and p n . 

• Compute the Kostant function K(G, h) for this system A + . 

• Compute a finite set F satisfying condition (C). 

• Compute a small element e such that p reg = p + e does not belong to any 
afline hyperplane of the form wX+H where H is a noncompact wall, w £ W c . 

• Compute for all w £ W c the number 



(11) 




wew c 



contribution. 



■IV 



N A +(WU - X- p n ) 



using the algorithm described in 13.6.11 
That is compute 

All w := V(w * (preg) — A,*4 + ) as explained in Figj6j 
Then compute 




G&F 



Me Allw 



ENDs 
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• Finally compute 



out := = e i w ) * contribution, 
tueWc 



Let us comment briefly: If w(u reg ) — A is not in the cone generated by non 
compact positive roots, the set All w is an empty set. In particular we may restrict 
the computation by diverse consideration to valid permutations which have 
some chance to give a non empty set, see Section I3~3l ) 

3.7.2 Symbolic 

The preceding calculation runs with symbolic parameter and we take advantage 
of this to find regions of polynomiality. 
Let's explain how. 

Let o be a chamber in t* for the system A of roots of g. Let U be the open set 
of (A, u) G a x t* such that wX — u does not belong to any non compact wall. Let 
(Ao,^o) £ U. Then we define R(Xq,uq) to be the closure of connected component 
of U containing (Ao,/^o)- This region is a cone in t* x t* with non empty interior 
and can be described by linear inequalities in A, u. 

For this domain we can compute a polynomial formula and state the following 
result. 

If A varies in a, the systems A+(A), A+(A) determined by A remains the same. 
We denote it by A+, A+. 

Furthermore, if (A, u) € R(Xq, /j,q), for any w G W c , the element wX — u lies in 
a tope t w for the system A+ which depends only of w. 

Theorem 42 The domain R(Xo, tq) is a domain of polynomiality for the Duistermaat- 
Heckman measure and thus is a domain of quasi-polynomiality for the multiplicity 
function m^. 

More precisely, for (A, u) £ R(Xo,uo) n (P r x i- > { r ), we have 
(12) mj= e(w)N T ™ t (wX-fi-p n ). 

The right hand side of Equation [12] is a quasi polynomial function of A, u, 
antisymmetric in u. It takes positive values if u is dominant for A^j~. Recall that, 
in this case, the multiplicity of u on A is the absolute value of the function 
above, (Sec l2.2p . 
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Of course, the symbolic calculation above, with the present approach, is limited 
to very small examples. 

Also, we are not able to determine the largest domains where the function m x 
is given by a quasi polynomial formula. 

3.7.3 Asymptotic directions 

We address now a simpler problem. We have the same setting that in the previous 
section, but we are now testing only the noncompact walls crossing in one fixed 
direction v. 

Let fio,Xo be given, with Ao G o n P g r an Harish-Chandra parameter, and 
Hq G a c n -P t r . Let v G a c be integral. We will do the calculation of when 
fH = Ho ~\~tv, with t > 0, is in the half- line in the direction v. In the application fio 
will be the lowest K-type Xq — p n of our discrete series tt x ° . We compute the values 
U where Ho + tv — wAo cross a non compact wall (other than the ones which may 
contain the line no + tv — w\q). These are the values where the line fit mav cross 
the domains of quasipolynomiality described above. We order this finite set of 
values < ti < t2 < ti < ■ ■ ■ < t s . Consider the interval Ii = [U, ii+i], < i < s, 
where to = and t s+ i = oo. 

Consider Ii n Z, an "interval" in Z, described by two integers [aj,6j], with 
Oj = ceil{ti) and 6j = floor{ti + \). 

The "interval" Ii n Z can also be reduced to a point. 

Then we find exponential polynomial function P l {t) on R such that m^ t is 
equal to P*(i) for t G J, nZ. 

If particular, u is an asymptotic direction, if and only if the last quasipolyno- 
mial P s (t) does not vanish. 

The algorithm is as follows. 

For each consecutive value ti,U + i, choose fir = fJ> + t r v with ti < t r < ti+i. 
Then, move very slightly fir i n H%- Then for each w G W c , w/i^. — Ao lies in a tope 
Tf for A+(A ). Then 

e HN^+ {A) (w/i t - A-p n ). 

The right hand side of this formula is an exponential polynomial function of 
t G Z. 

The algorithm implementing this procedure is described in FigJHJ Let us remark 
that our algorithm implementation is for type A r and thus the Pj are polynomials. 
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4 Blattner's formula for U(p,q) 



4.1 Non compact positive roots 

With the notation of Section [JJ we let G = U(p, q) and K = U p x U q be a maximal 
compact subgroup. 

Let E be a p + ^-dimensional vector space with basis (i = 1, . . . , p + g) and 
V as in Ex. [To]) , Let r = p + q — 1. Consider the set of roots 

A = ±{ej — ej 1 1 < i < j < p + q}. 

We then choose T to be the diagonal subgroup of U(p,q), and identify t* with 
E. In this identification the lattice of weights is identified with W +q : the element 
(m, . . . , n p+q ) giving rise to the character 

t = (exp(i0i), . . . ,exp(i0 p+(z )) -> e inidl ■ ■ ■ e in "+" e "+^. 

The system of compact roots A c is 

A c = ±{ei - ej 1 1 < i < j < p} U ±{ei - ej \ p + 1 < i < j < p + q}. 
The system of non compact roots is 

A n = ±{ei - ej | 1 < i < p,p + 1 < j < p + q}. 

Let A be the Harish Chandra parameter of a discrete series for G and \x the 
Harish-Chandra parameter of a finite dimensional irreducible representation of K. 

Because discrete series are equivalent under the action of the Weyl group of K, 
then we may assume that A = [a, 0\ where a = Y^=i a i e i = [ a i> • • • > Q p]> a i > 
a 2 > • • • > a p and (3 = Yn=p+i Pi e i = \Pl> ■ ■ ■ > Pq]i Pi > fa ' ' ' > Pq- 

Here cq, /3j are integers if p + q is odd, or half-integers if p + q is even, that is 
we fix as system of positive compact roots the system A+ = {e, — ej \ 1 < i < j < 
p} U {ej - ej \ p + 1 < i < j < p + q}. 

We parametrize fi G P t r by another couple 

H := [a,b] = [[01,02,..., Op], [&!,...,&,]] 

with a% > ■ ■ ■ > a p and b± > ■ ■ ■ > b q . 

Here are integers if p is odd, half-integers if p is even. Similarly bj are 
integers if q is odd, half-integers if q is even. 

As the center of G acts by a scalar in an irreducible representation, we need 
that the sum of the coefficients of A has to be equal to the sum of the coefficients 
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of ji for the multiplicity of fi in ir x to be non zero. Thus A — \i is in V, (see Ex. 
EU). 

We now parametrize the different dominant chambers of t* modulo the Weyl 
group of K by a subset A of [1, 2, . . . , r + 1] of cardinal p. Let B its complementary 
subset in [1, 2, . . . , r + 1]. 

To visualize A, B we write a sequence of lenght p + g of elements a, 6 with a 
in the places of A, b in the places of B: for example if A = [3, 5] and B = [1, 2, 4], 
then we write [b, b, a, b, a] or simply bbaba. Now we use this visual aid and describe 
a permutation wa of the index [1, 2, . . . ,p + q], by putting the index [1, . . . ,p] in 
order and in the places marked by a, and the remaining indices [p + 1, . . . ,p + q] in 
order and in the places marked by b, precisely wa '■ [1, 2, 3, 4, 5] — > [3, 4, 1, 5, 2]. The 
elements wa where A varies describe a system of representatives of S p+(? / (S p x £ g ), 
(£„ being the permutations on n letters), that is also the chambers of t* for A(g, t) 
modulo W c . In the above the chamber is described by {h = [hi, h 2 ,h^,h/i, h§\ h% > 
hi > hi > h$ > /12} 

Let a s tandard be the chamber ai > a 2 > • • • > a p > (5i > (3% • • • > fiq- 

Then if A G wa& standard, w e have A+(A) = A+ and A+(A) is isomorphic 
to A + (A, B), by relabeling the roots via w^ 1 ■ The next example will clarify the 
situation. The subset A can be read from A: we reorder completely the sequence 
A and define A as the indices where the first p elements of A are relocated. 

Example 43 Let G = £7(2,3) with compact roots A c = ±{ei — e2,e3 — e^,e^ — 
e§, e^ — e^} and noncompact roots A n = ±{ei — e^, ei — e^, ei — e^, e<i — e2 — e±, e2 — 
65}. Let A = [a, f3] with a = [4, 2] = 4ei + 2e2 and (3 = [6, 5, 3] = 6es + 5e4 + 3e5. 

Then A = [3,5], B = [1,2,4] that is the configuration bbaba. The system of non 
compact positive roots for A is — e\, e± — e — 1, e — 1 — es, e^ — e<i, e^ — e2, e§ — e^, 
isomorphic to A + (A, B) by the relabeling of the roots suggested by w^ 1 , that is 
e3 = /1, e 4 = f 2 , ei = / 3 , e 5 = / 4 , e 2 = /s- 

Thus relabeling the roots, our calculations will be done for A + (A,B) inside Af~, 
where A + (A,B) is given in Example 1161 

Remark here that A + (A, B) is irreducible if p and q are strictly greater than 1. 
In contrast, when p or q = 1, the system is fully reducible. Consider for example 
the case p = 1. 

Example 44 

In this case A has only 1 element and the system A + (A, B) has r elements and 
is a base of V. Thus A + (j4, B) is fully reducible in the direct sum of p + q — 1 one 
dimensional systems. 

For example take £7(1, r) with A = [1] and B = [2, . . . , r+1]. Then A+(A, B) = 
{ei — eg, ei — e%, . . . , e\ — e r+ {\ is isomorphic to Af x A^ x • • • x . 
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Remark that when A±, A2 have the same number of elements, although the sys- 
tem of noncompact roots ±A^{Ai, B\) and ±A+(A2, B2) are clearly isomorphic, 
the combinatorial properties of A+(A, B) may vary. 

For example, (see Figure [5]), if A = [1,2], B = [3,4], the cone generated by the 
non compact roots has basis a square and is not a simplicial cone. If A = [1,3] 
and B = [2, 4], then the cone generated by the non compact roots is the simplicial 
cone generated by e± — e2 — e%, — e±. 

4.2 Algorithm to compute MPNS: the case of A+(^4, B) 

With the notations of Exfl6| we denote by A a proper subset of [1, 2, . . . , r + 1] 
(with r = p + q — 1) and by B the complementary subset to A in [1, 2, ... r + 1]. 

Given v € V, and not on any admissible hyperplane, we describe the algorithm 
to compute V(v, A+(A, B)). 

If p or q = 1, roots a in A + (^4, B) form a basis on V, thus there is only one 
maximal nested set M = {{a}, a € A + (A, B)}. Thus V(v, A + (A, B)) is empty or 
equal to {M} depending if v belongs to the cone generated by A + (A, B), or not. 
This is very easy to check. 

If p > 1 and q > 1 , we determine the set V(v, A + (A, B)) by induction, going 
to admissible hyperplanes. 

If L C [1, 2, . . . , r + 1] is a proper subset of [1, 2, . . . , r + 1], we will also use 
the notation V = [i ^ L\l < i < r + 1] for the complement of L. We denote by 
Hl ■= {v £ V I YlieL v i = 0} the hyperplane determined by L; the hyperplane 
Hl is equal to the hyperplane Hl 1 determined by V . 

It is very simple to describe A + (A, l?)-admissible hyperplanes, that is noncompact 
walls. The description is an adaptation of the A+-admissible hyperplanes that 
appear in [3]. 

Keeping A fixed, with \A\ ^ 1, r, we consider hyperplanes Hl indexed by 
subsets L C 1, 2, . . . , r + 1] with the following properties: 

• if \L\ 7^ 1 or r, then i^x, is a noncompact wall if and only if both A and B 
intersect L and L'. In this case A+(A, B)(~]Hl is the product of two systems 
A+(A n L, B n L) x A+(yl n L' , B n L') and thus reducible. 

• if L is of cardinal 1, then if/, is a noncompact wall. In this case A+(A, B) n 
i?L is A+(AnL',Bfl L') and thus irreducible. 

At this point to compute the MNPS or better, as we explained the M's, we 
can proceed as in Figj6j The algorithm is outlined in Figd 

We conclude with the following observation. A necessary and sufficient con- 
dition for the set MPNS(v, A+(A, B)) to be non empty is that v belongs to the 
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cone generated by A+(A,B). As far as we know, the equations of this cone are 
not known, except in a few cases. It is clearly necessary that v belongs to the 
simplicial cone generated by all positive roots. To speed up the calculations, we 
check this condition at each step of the algorithm. 

We conclude with a simple example with p = 2,q = 2 and A = [1,2], B = 
[3,4]. We follow the outline described in FigJSJ The highest non compact root 
is 6 := ei — e^. There are 3 noncompact walls not containing the highest root. 
L = [1], [4], [1, 3] We choose a vector v = [4,3,-2,-5] not on any noncompact 
walls. Then [1], [4], [1, 3] are all such that v and 9 are on the same side. 

For L = [1], the v projection do not belong to the cone generated by A^(A, B)(~) 
Hl = [&2 ~ e3,e 3 - e 4 ]. 

For L = [4], we obtain the element M := {[1,2,3,4], [1,3], [2,3]} in V(v, A+(A, B)). 

For L := [1,3], we obtain the element M = {[1,2,3,4], [1,3], [2,3]} in V(v, A+(A, B)). 

4.3 Valid permutations 

Let w S W c . Remark that if wfi — A does not belong to the cone of non compact 
positive roots, then the term corresponding to w in Blattner formula is equal to 
0. It is important to minimize the number of terms in Blattner formula. To this 
purpose, we use a weaker condition: we say that w € W c is a valid element if 
wfj, — A is in the cone spanned by (all) positive roots. Thus if w is not valid, the 
corresponding term to w in Blattner formula is equal to 0. As there is a simple 
description of the faces of cone spanned by all positive roots (it is the simplicial 
cone dual to the simplicial cone generated by fundamental weights), there is a 
simple algorithm that constructs valid permutations one at the time depending on 
the conditions they have to satisfy, instead of listing all the elements of W c . The 
corresponding algorithm is used in [3|,[B], and we just reproduced it. 

5 Examples 

Example 45 

We consider the discrete series representation indexed by A and we test for the 
multiplicity where \i is a K type. We write ^i owes t for the lowest if -type. We 
use the algorithm whose command is : 

>discretemult(A,/i,p,q) 
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m^: numeric case 


Group 


Input 


Output 


Time 


U(3,3) 


A=[[31/2,15/2,9/2], [5/2,3/2,-5/2]] 






*W«t=U17, 9, 6j, [1, 0, -4]] 


1 


0.026 sec. 


/i=[[1017, 1009, 1006], [-999, -1000, -1004]] 


9 


0.97 sec. 


M =[[100017, 10009, 10006], [-9999, -10000, -100004]] 


9 


0.91 sec. 


U(3,3) 


A=[[31/2, 19/2, 11/2], [ 15/2, 7/2, -37/2]] 






No W est = m, 11, 6], [7, 2, -20]] 


1 


0.073 sec. 


M= 


[1017, 1011, 1006], [-993, -998, -1020JJ 


275 


0.529 sec. 


M= 


100017, 10011, 10006], [-9993, -9998, -100020JJ 


11700255 


0.538 sec. 


U(4,4) 


A= 


11/2,7/2,3/2,-1/2], [9/2,5/2, 1/2,-3/2]] 






rt i Ues i = [[15/2, 9/2, 3/2, -3/2], [11/2, 5/2, -1/2, -7/2JJ 


1 


0.565 sec. 


^=[[2015/2, 9/2, 3/2, -3/2], [11/2, 5/2, -1/2, -2007/2JJ 


120495492015 


3.493 sec. 


U(4,4) 


A=[[ll/2, 9/2, 7/2, 5/2], [3/2, 1/2, -1/2, -3/2]] 






rto» esi = [[15/2, 13/2, 11/2, 9/2], [-1/2, -3/2, -5/2, -7/2]] 


1 


0.334 sec. 


/u=[[20015/2, 2013/2, 211/2, 29/2], [-21/2, -203/2, -2005/2, -20007/2JJ 


1 


273.719 sec. 


U(5,4) 


A=[[5, 3, 1, -1, -3], [4, 2, 0, -2]] 






M io » esi = [[7, 4, 1, -2, -5], [11/2, 5/2, -1/2, -7/2]] 


1 


3.952 sec. 


M =[[1007, 4, 1, -2, -5], [11/2, 5/2, -1/2, -2007/2]] 


120495492015 


13.752 sec. 


U(5,5) 


A=[[ll/2 J/2,3/2,- 1/2,-5/2], [9/2,5/2,1/2,-3/2,-7/2]] 






AW s t=[[8, 5, 2, -1, -4], [6, 3, 0, -3, -6JJ 


1 


51.910 sec. 


/u=[[106,4,2,0,-102], [104,2,0,-2,-104]] 


1458704380546472381 


163.104 sec. 



Example 46 

We consider the discrete series representation indexed by A, a direction v and we 
test for the multiplicity rn^ +t ^ where /x = Slowest is the lowest -fT-type. We use the 
algorithm whose command is 

> function discrete mu direction lowest- ( A,w,p,q) 

For completeness we list ^lowest relative to each example. 
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m u+tv i * G asymptotic case 


Group 


Input 




Output 


Time 


U(2,3) 


A=[[9, 7], [-1, -2, -13]] 








n lowest =[[21/2, 17/2J, [-2,-3, -14JJ 


tT=[[l, 0], [-1, 0, 0]] 


ra A . = - 


1 if t = 
if t > 1 


Al 


0.054 sec. 


?=[[6, 1], [-1,-1,-5]] 


m A , , - = < 


1 if t < 10 
if t > 11 


A2 


0.736 sec. 


v=[[l, 0], [0, 0, -1]] 


m A , = 1 if t > 


A3 


0.26 sec. 


U(2,3) 


A=[[59, 39], [51, 7, -156]] 








W «,e«t=lll21/2, 79/2], [51, 6, -157JJ 


«=[[1,0], [0,0,-1]] 


m A ,,^ = t + 1 


A4 


0.71 sec. 


U(2,4) 


A=[[341/2, 49/2], [-3/2, -5/2, -11/2, -371/2]] 








rto» es t = [[345/2, 53/2], [-5/2, -7/2, -13/2, -373/2]] 


tT=[[6, 1], [-1,-1, -1, -4]] 


A 1 1 if t < 2 

m M +*« ~ \ if t > 3 


A5 


46.754 sec. 


U(3,3) 


A=[[343/2, 31/2, 21/2], [-13/2,-19/2,-363/2]] 








rt otues t = [[173, 17, 12], [-8,-11,-183]] 


v=[[6, 1, 0], [-1, -1, -5]] 




' 1 if t = 
3 if t = 1 
6 if 2 < i < 171 
3 if t = 172 
1 if t = 173 

. if t > 174 


A6 


52.020 sec. 


v=[[l, 1, 0], [0,-1,-1]] 


\ I t + 1 if < i < 155 
m+*s \ 156 if t > 156 


A7 


7.1730 sec. 



where 



Al 


= [[-(1/2) * + (1/2) * t + 1, [0, Ojj, [1 + (1/2) * - (3/2) * t, [1, ljj, [0, [2, in/]]] 


A2 


= [[-(15/2) * ^ + (7/2) * t + 1, [0, 0]], [1, [1, 10]], [66 + (1/2) * - (23/2) * t, [11, 11]], [0, [12, inf]]] 


A3 


= [[-(1/2) * + (1/2) * t + 1, [0, OJJ, [1, [1, infl\\ 


AA 
A5 


= [[t + l,[0,m/]] 

= [1 + (10/3) *t- (7/2) * t 2 + (25/6) * t 3 , [0,0]], [1 + (1/3) * t - (1/2) * t 2 + (1/6) * i 3 , [1, 1]], [1, [1,2]], 
[6 - (7/2) * t + (1/2) * t 2 , [2, 3]] , [0, [3, inf]] 


A6 


= [1 + (13/4) *t- (51/8) * t a + (61/4) * t a - (89/8) * t 4 , [0,0]], [1 + (3/2) * t + (1/2) * t 2 , [1, 1]], 
[-3 + (27/4) * t - (1/8) * t 2 - (3/4) * t 3 + (1/8) * t 4 , [2,3]], [6, [4, 170]], 
[32664996 - (9380059/12) * t + (168227/24) * t 2 - (335/12) * t 3 + (1/24) * t 4 , [171, 172]], 
[15225 - (349/2) * t + (1/2) * t 2 , [172, 173]] 

[78155000 - (10718575/6) * t + (183749/12) * t 2 - (175/3) * t 3 + (1/12) * t 4 , [174, 175]], [0, [176, inf]] 


Al 


= [1 + (5/4) *t+ (1/24) * t 2 - (1/4) * t J - (1/24) * t 4 , [0,0]], [t + 1, [1,154]], 

[23726781 - (2457885/4) * t + (143207/24) * t 2 - (103/4) * t 3 + (1/24) * t 4 , [155, 156]], [156, [156, inf]] 



Remark that in some of the examples above, it can happen that although the 
polynomials Pi and Pj+i (displayed in the last table giving the Ai) are different, 
the polynomial Pi+i may coincide with Pi on ij (recall that several polynomials 
can have the same values on I,nN). Thus in this case, we join the two intervals 
Ii and /j+i and give only the polynomial Pj+i. This streamlining of the function 
m ^,+tv ^ s gi ven m the third column of the table describing the asymptotic behavior 

° f mX a+tv- 

We now give an example on [7(3, 4). 
Example 47 
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discrete:=[[473, 39, 1], [3, 51, 5, -572]]; 
directions, [[1, 0, 0], [0, 0, 0, -1]]; 



>f unction_discrete_mul_direction_lowest (discrete , direct ion, 3,4) ; 

[1+51/20 t-1/120 t-5-1/360 t'6+851/360 t'2+23/24 t~3+5/36 t'4,[0,0]], 
[1+31/12 t+19/8 t'2+11/12 t"3+l/8 t"4,[l,37]], 

[-3262622+2687514/5 t+73/240 t'5-1/720 t"6-13275857/360 t'2+64795/48 t~3-3977/144 t~4, [38,39]] , 
[-265030+27790 t-1090 t'2+20 t"3, [39,44]] , 

[-9631849+79305707/60 t+29/80 t'5-1/720 t"6-3399664/45 t'2+110609/48 t~3-5675/144 t~4, [45 ,45] ] , 
[27182687-212385511/60 t-31/40 t~5+l/360 t'6+69073219/360 t'2-132929/24 t~3+1619/18 t~4, [46 ,46] ] , 
[-784945+886169/12 t-20959/8 t'2+511/12 t'3-1/8 t~4, [47,83]] , 

[469370132-337238937/10 t-167/240 t"5+l/720 t'6+363465857/360 t"2-774005/48 t'3+20897/144 t"4, [84,85]] , 
[5790400-235000 t+2820 t"2 , [86 , inf ] ] 

Thus the multiplicity m^ +t ^ can be completely described by the following 
piecewise polynomial function: 

(1 + (31/12) * t+ (1/8) * t 4 + (11/12) * t 3 + (19/8) * i 2 if < t < 39 

-265030 + 27790 * t + 20 * t 3 - 1090 * t 2 if 40 < t < 46 

-784945 + (886169/12) * t - (1/8) * t 4 + (511/12) * t 3 - (20959/8) * t 2 if 47 < t < 85 

5790400 - 235000 * t + 2820 * t 2 if 86 < t. 

The time to compute the example is TT := 19.487 and the formula says 
for instance that, for A =discrete and \i = [[475, 40, 0], [103/2, 9/2, 5/2, -1147/2]] 
the lowest K-type, then 

TO^nnnnnnrw = 1 127995300005790400. 



6 The program: "Discrete series and K mul- 
tiplicities for type A r " 

We give a brief sketch of the main steps for the algorithms involved in Blat- 
tner's formula. 



6.1 MNPS non compact 

We outline the algorithm that computes directly M for M G V(v, A + (A, B)). 
We are taking advantage of the fact that we know the M's in the case of 
\A\ = 1, r, as we saw in Ex. HH In the following scheme p, q are integers, A C 
[1, 2, . . . , p + q] is a set of cardinality p, B is the complement subset defining 
U(A, B), 8j is the highest noncompact root for /. If L C [1, 2, . . . ,p + q] we 
denote by V the complement set. 
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Input [v, A, I] , v a vector and A C / = [1, 2, . . . , p + q], \A\ = p 
proceed by induction on the cardinality of A. 

if |A| = 1 or \A\ = r write the unique M, M E MNPS determined by the situation 

if v 6 C(M) then the output is M, 
construct the hyperplane Hl- 
check if v and 9j are on the same side then Hl 

if not skip the hyperplane 
define the projection v' = proj^ (v) of v on Hl along 8j 
compute [v\ , Ai , Ji], [V2, Aq., h] 

where A\ = L n A, A2 = L' n A, h = L, h = L' and 

vi , V2 are the components of v' on L and L' respectively. 

if t>i(resp. V2) is not in the positive cone for A(7i) then skip the hyperplane 
if = 1, apply the induction and compute Mi, Mi e MNPS(vi , Ai , h), 
add to Mi the root dj x (do the same if | ^4.2 1 = 1 ) 
else apply the induction and compute Mi, Mi £ V(vi, Ai, Ii),i = 1,2 
do the cartesian product M\ X M2 and add to each set the root 81 
collect all M's, for the wall L 

end of loop running across L's 

end induction 

return the set of all M's, M £ MPNS for all hyperplancs 



Figure 7: V(v,A+(A,B)) 



6.2 Numeric 

The scheme is described in FigJHl 

6.3 Asymptotic directions 

We fix the parameter A and /i , regular in the chambers a, a c and a weight v. 
We want to compute rn^° +tiI . In the application /i will De the lowest K-type. 
The scheme is described in FigJHl 
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Subroutines: 

• Procedure to find .4+ -admissible hyperplanes. 

• Procedure to deform a vector: DefVec n c(v, A^) 

• Procedure to compute M, M £ V(v, A + ) as in FigU 

• Procedure to compute Kostant function K(h) = K(0,h) = -j=j e _ or more generally 

K(g,h). 

• Compute the valid permutation Valid(u, v) C W c 



Input: A, fi 

Compute A + = A^(A): 
Compute Valid(X, fi) 

for each w £ Valid(X, fj,), 

compute fireg = DefVec nc (v,A + ) 

compute All w := r P(wfj, reg — A, A J) 

compute cont w = J2jft cAll ^ les j^ K ( w ^ - A - p n ) 

end of loop running across w's 
collect all the terms and return 

™» = E m6 y a iid(A, M ) <w)cont w 

Figure 8: Blattner's algorithm (numeric case) 

Subroutines: 

• Procedure to find A + -admissible hyperplanes. 

• Procedure to deform a vector: DefVec nc (v,A + ) 

• Procedure to compute M, M 6 \CP(v, A + ) as in Fig[7] 

• Procedure to compute Kostant function K(h) = K(0,h) = y=j e ( 1 _ e -e») or more generally 

K{g,h). 



Input Ao and v for each H noncompact wall 
if (H, v) = then skip H 

else if (H, Ao — wfio)(H, v) < then skip H else 
collect tii = (H, Ao — wfiQ)/(H, v) 
end of loop running across H's 

order list t' H s as [to, ti . . ,t 3 ] where to = 0, t s = oo 
choose an interior point ti in each interval [tj,i i+i] 

Compute polynomial on each [t^ , ii+ i ] following the scheme FigOand the ordered basis determined by ti 
output: the sequence of values ra^° +tg valid on [tj,tj_|_i], tgW 



Figure 9: Blattner's algorithm: asymptotic case 
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